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COMPUTER  AIDED  APPROACH  TO  THE  DESIGN  OF  Y-JUNCTION 
STRIPLINE  AND  MICROSTRIP  FERRITE  CIRCULATORS 


1.  INTRODUCTION 

New  interest  is  evident  in  modernization  of  ferrite  component  technology.  Monolithic 
microwave  circuit  technology,  which  was  perfected  in  the  1980s,  has  permitted  the 
development  of  multidecade  bandwidth  amplifiers,  multioctave  bandwidth  subsystem 
components,  and  full  systems  that  have  more  than  an  octave  of  bandwidth.  The  Y-junction 
circulator  is  the  workhorse  of  ferrite  components.  It  can  be  used  as  a  3-port  circulator,  2- 
port  isolator,  or  a  single-pole,  double-throw  switch.  For  new  circulator  technology  it 
would  be  useful  to  increase  the  bandwidth  capability  (currently  about  an  octave),  minimize 
the  size  (perhaps  by  using  matching  circuits  other  than  quarter- wavelength  transformer  s), 
and  create  design  approaches  that  are  compatible  with  monolithic  circuits.  The  first  step  is 
to  have  on  hand  an  accurate  theoretical  analysis  tool  to  calculate  intrinsic  circulator 
performance,  thereby  permitting  simultaneous  circulator  and  matching  circuit  optimization. 
The  following  is  an  effort  to  assemble  a  concise  circulator  compendium  that  consists  of 
background  and  tutorial  information  on  ferrite  media,  wave  propagation  in  ferrite  media, 
junction  circulator  theory  and  evolution,  modem  improvements  in  the  theory  to  permit 
wideband  analysis,  and  finally  a  computer  program  to  perform  the  wideband  analysis. 

The  landmark  theoretical  works  on  Y-junction  circulators  are  those  of: 

•  Bosma  [1]  (solves  the  simplified  boundary- value  problem  of  the  circulator  ferrite 
disk  by  using  a  Green’s  function  approach  that  relates  the  axial  component  of  the 
electric  field  to  the  circumferential  component  of  the  magnetic  field  at  the  perimeter 
of  the  disk); 

•  Fay  and  Comstock  [2J  (describe  the  operation  in  terms  of  counter-rotating 
propagation  modes  in  the  ferrite  disk); 

•  Wu  and  Rosenbaum  [3]  (predict  octave  bandwidth  microstrip  circulator  operation 
theoretically  using  Bosma's  Green's  function  analysis);  and 

•  Schloemann  and  Blight  [4]  (rederive  the  Bosma  form  of  the  Green's  function  for 
\ieff  <  0).  When  the  bias  field  is  low. 


Veff=  (V-2  -  *2)/  H 


(1) 


is  real,  where  p  and  k  are  the  off-diagonal  components  of  the  ferrite  material's 
permeability  tensor.  Mu-effective  (\ieff)  becomes  real  and  negative  at  frequencies 
well  below  the  circulator  design  center  frequency,  at  and  below  the  point  where 


/=  (y/2ic)  •  4jc  Ms, 


(2) 
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in  which/is  the  frequency,  y  is  the  gyromagnetic  ratio  of  electron  magnetic 
moment  to  angular  momentum,  (yUn)  is  approximately  2.8  MHz/Oersted,  and 
4 nMs  is  the  ferrite  material's  saturation  magnetization. 

The  new  theoretical  work  reported  here  includes  the  derivation  of  the  analytical  link 
between  the  results  of  Refs.  [3,4],  the  inclusion  of  the  imaginary  part  of  the  dielectric 
permittivity  of  the  ferrite  material  in  the  calculations,  and  the  development  of  the  computer 
program  for  quantitative  analyses.  The  first  item  facilitates  the  use  of  a  set  of  published 
equations  for  part  of  the  computer  programming;  and  the  second  permits  the  effects  of 
dielectric  losses  in  the  ferrite  material  to  be  calculated. 

Some  of  the  most  important  experimental  stripline-junction  circulator  work  has  been 
that  of  Simon  [5],  which  gives  the  results  of  tests  on  a  wide  range  of  materials;  Salay  and 
Peppiatt  [6],  which  clarifies  the  sign  of  the  reactive  part  of  the  input  impedance;  and 
Schloemann  and  Blight  [4],  which  shows  very  wideband  performance  over  3  :1  bands. 

2.  THEORY 

(a)  Ferrite  Medium 


A  biased  ferrite  material  has  complex  permittivity  £  and  tensor  permeability  [l,  and  its 
elements  may  be  complex.  The  tensor  character  of  the  permeability  is  developed  from  the 
Landau-Lifshitz  equation  of  motion  of  the  magnetization  vector  M  in  the  material 


dM 

a  t 


-  y  (MxH)  + 


ya  (MxH)  x  M 


m 


(3) 


with  a  small  RF  signal  applied.  M  is  magnetization,  H  is  magnetic  field,  and  a  is  an 
experimentally  determined  parameter  that  adjusts  the  magnitude  of  the  bracketed  damping 
term.  In  rectangular  coordinates,  with  the  magnetic  bias  field  in  the  z-direction,  H  and  M 
are 


H  =  ihx  +  jhy  +  iiHB 


(4) 


and 


M  =  i  mx  +  j  my  +  k  Afjg  , 


(5) 


where  mx  «  Mb  and  my  «  Mb  ,  with  the  lower  case  letters  indicating  the  RF  values  and 
the  capital  letters  indicating  the  bias  values.  (Usually  Hg  is  large  enough  to  make  Mg  equal 
to  the  saturation  value  M$.)  The  solution  has  the  form 


j^MbHb  ^(o+XHbIMb) 

h  - 

y[  Mb  {j(o  +  XHb  /  Mb)  -  X  Hb] 

HB+[j(d  +  XHb  !  Mb f" 

nXQ 

y2  Hg  +{j(o  +  XHb  /  Msf 

hyQ,  (6) 


where  A  =  -y  a  I M  |  and  I M I «  Mb,  with  the  zero  subscripts  indicating  that  the  harmonic 

time  dependence  is  not  included  in  the  values  (wi*  =  By  convention,  Eq.  (6)  is 

written: 

flijto  =  Xxx  hXQ  -j%yx  hy0  ,  (7) 

and  the  similar  solution  for  myQ  is : 
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myo  —  J  Xxy  hx o  +  Xyy  • 


In  the  solutions  Xyy  =  Xxx  and  Xyx  =  -Xxy  so  that  the  intrinsic  RF  magnetic  susceptibility 
tensor  can  be  written  as 

Xxx  ~JXxy  0 

X=  jXxy  Xxx  0  ,  (9) 

0  0  0 

*“  mJ 

assuming  mz  is  negligibly  small  relative  to  Ms-  In  the  CGS  electromagnetic  system  of 
units  used  in  this  pan  of  the  derivation, 


where 


b=h  +  47tm  =  (l+4^x)h  =  fih. 


'l  0  op* 

h=  0  1  0  hy 

.0  o  iJa 


1  +  4  X**  -fix  Xxy  0 

4  =  fix  Xxy  \  +  An  Xxx  o  ,  (12) 

.0  0  1 

which  may  be  written  as 

[  *1  -jK  0  ' 

4  =  4  0  .  (13) 

.0  0  1 

This  is  the  intrinsic  permeability  tensor  of  the  ferrite  material  and  is  related  to  material 
parameters,  bias  field,  and  frequency  as  shown  in  Eq.  (6), 


where 


4=1+4  7tXxx=tf-jtf', 


K  =  4  x  Xxy  =  K'-jK". 

(b)  Plane  Waves  in  Magnetized  Ferrite  Media 

For  an  RF  field  in  any  linear,  isotropic,  homogeneous  medium 

VxE=-^  | 
dt  l 

Vxh  =  -^- 
dt 


with  D  =  €  E  and  b  =  £  h  .  Consider  plane  waves  of  the  form 
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(E  ,  h  ,  b  ,  D  )  =  (  Eo  ,  h0  ,  bo  ,  Do  )  e/  *> ' -  r*»  •  r  ,  (17) 

where  the  0  subscripted  values  do  not  have  time  and  coordinate  dependence,  where  n  is  a 
unit  vector  in  the  direction  of  plane  wave  propagation,  and  r  is  related  to  the  system 

coordinates  (  r  =  ix  +  jy  +  kz  in  rectangular  coordinates).  Solving  Eq.  (16)  in  magnetic 
«-» 

terms  that  involve  It  yields 


a2  e  b0  =  T2  n  x  (  n  x  h0  )  (18) 

which  may  be  solved  for  T,  the  propagation  constant. 

The  propagation  constant  is 

r±  =jarfepo 


l (pi2  -pi  -  k2)  sin2  0  +  2 pi± [(pi2  -pi  -  k2)  sin*  0  +  4  k2  cos2  o]1/2 
\  2  [0z-l)sin2  0  +  l] 

The  term  inside  the  braces  is  the  relative  effective  permeability  that  depends  on  0,  the  angle 
between  the  direction  of  propagation  and  the  magnetic  bias  field,  as  well  as  on  the  material 
parameters  contained  in  £  =  £1)  £>  ,  piQ  ,  pi ,  and  K,  and  on  frequency  contained  in 
1 0  =  2  icf .  For  propagation  perpendicular  to  Hg,  as  is  the  case  in  a  junction  circulator 
"puck"  of  ferrite  material, 


\l/2 


2/1 


=ja)SefI0{pieff}112 .  (20) 


Then 


=  K2)/pi,  (21) 

which  corresponds  to  a  linearly  polarized  wave  with  the  h  field  perpendicular  to  Hg.  Also, 

Veft-  =  1.  (22) 

which  corresponds  to  a  linearly  polarized  wave  with  the  h  field  parallel  to  Hg.  It 
propagates  as  if  in  air  except  for  permittivity  effects.  In  the  junction  circulator,  pi^ff+  is  the 
operative  term. 

The  preceding  theoretical  discussion  is  sufficient  background  to  proceed  to  the  special 
problem  of  Y-junction  circulators.  General  information  on  additional  theoretical  matters  and 
on  other  kinds  of  ferrite  applications  is  provided  by  Soohoo  [7],  and  by  Linkhart  [8]. 

(c)  Stripline  and  Microstrip  Junction  Circulators 

The  field  problem  for  the  3-port  symmetrical  Y-junction  stripline  circulator  has  been 
addressed  by  Bosma  [1]  by  using  the  notation  shown  in  Fig.  1.  The  basic  assumptions  in 
the  analysis  were: 
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+=n 

ISOLATED 


INPUT 

Fig.  1  -  Circulator  coordinate  reference 

(1)  Stripline  configuration:  Two  equal  thickness  ferrite  disks  between  two  ground 
planes  and  a  zero  thickness  center  conductor  disk  of  the  same  radius  with  three  equally 
spaced  coupling  lines  of  width  w,  as  shown  in  Fig.  1.  Bosma  notes  that,  as  is  the  case  in 
striplines,  the  electromagnetic  fields  are  equal  in  magnitude  and  opposite  in  polarity  on 
opposite  sides  of  the  center  conductor,  so  that  the  field  problem  in  the  ferrite  need  be  solved 
for  only  one  disk.  Therefore,  his  solution  is  approximately  valid  for  microstrip. 

(2)  The  fields  in  the  disk  do  not  depend  on  the  z  coordinate,  so  that  the  solution  need 
be  only  2-dimensional  (2-D). 

(3)  Except  at  the  connections  of  the  coupling  striplines,  no  radial  current  can  flow  from 
the  edge  of  the  center  conductor. 

(4)  The  tangential  component  of  H  <{,  (R  ,  <)>)  in  the  ferrite  is  equal  to  the  magnetic  field 
intensity  in  the  corresponding  coupling  stripline  at  the  same  place. 

(5)  H  (R  ,  40  is  constant  over  the  coupling  stripline  width  w,  and  it  is  0  elsewhere. 

In  the  cylindrical  coordinate  system,  the  values  of  H<|,  (R  ,  <j>)  are  Ha  at  the  input  port,  Mb 
at  the  output  port,  Hc  at  the  isolated  port,  and  0  elsewhere.  In  general  Ha,  Hb ,  and  Hc  are 
complex. 

(6)  The  electric  field  intensity  in  the  disk  has  a  z  component  only,  Ez  (r,  <j>),  or 
EZ(R ,  <j>)  at  the  disk  perimeter. 

In  Ref.  [1]  Bosma  indicates  that  Ez  (r,  <)>)  satisfies  the  Helmholtz  vector  wave  equation 

V2E  =/i€  —  (23) 

dt2 


or  assuming  harmonic  time  dependence 
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V2E  =  -(o2ne  E  . 


(24) 


So,  in  cylindrical  coordinates 


a2#,  i  dE,  i  c^Er  ,  _ _ ,2  ^ 

^+^+^7J+i<0'rm  £’ 


=  0. 


From  Maxwell's  equation  assuming  harmonic  time  dependence, 

VxE  =  jcofMM  H 


comes 


giving 


and 


idEz 

r  a<{> 
be2 

dr 

0 


=jCW 


M-  ~j<  0 

JK  n  0 

0  0  1 


Hr  = 


■1  dEz  K  dEz 
■^a< >  n  dr 
WoHeff 


//$  = 


ar  ^  a<> 


'  Hr  ‘ 

//$ 

. 

Hz 

(25) 


(26) 


(27) 


(28) 


WoV-eff 

A  Green's  function  is  then  introduced  that  relates  //<j>  and  Ez  at  the  edge  of  the  disk: 


(29) 


,(*,<►)  =  [  G(R$  ;  R,V)  H^Rtf)  d*'. 
J-n 


(30) 


For  <j>  in  radians, 
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Ha  at(-|-T)<(|.<(-|  +  'P) 

Hb  at(|-4<)  <$<  + 

Hc  at(jr- «§')  <<><  (I  +  4-) 

Zero  everywhere  else 

Then,  using  G  (<j>  ;  <J>’)  for  G (R  ,<j)  ;  /?,< {>’),  and  for  small  \|f , 

H*  ,  - f )  =  £,  =  2  4-  [C(-  f ;  -  f  )».  ♦  C(-  * ;  f  )K*  +  c(-  X ;  ,)HC] 

f)=E«=  2  4<[c(|;-|)H0  +  C(|;f  )Hi,  +  C(|;)t)ffc]  .  (32) 

and 

EJtR,n)  =  Ec  =  2'¥  [g(jc  ;  -  |)ffa  +  G(k  ;  |)tf +  G(k  ;  7t)//c]  J 

Note  that  £4,  Eb,  and  Eq  are  the  electric  field  intensities  at  each  of  the  three  circulator 
ports:  input,  output,  and  isolated,  respectively. 

Bosnia  [1]  notes  that  a  solution  to  Eq.  (25)  is  a  series  with 

Ez,n(rA)  =  anJn(kr)ejn* ,  (33) 

where  Jn  is  the  nth  order  Bessel  function  and  k  is  the  propagation  constant 
cofjiQ  Eft  Heff  £f) where  //<■#■  and  prefer  to  the  ferrite  material.  Then  from  Eq.  (29), 

#*>,<!>)  =  j  ^ '  (34) 


Finally,  the  Green's  function  is 
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G(M>  ;/?,<{>•)  =X  +  Y  ) 

X  _  ~j  Zeff  Jpikr) 

2  K  JQ(kr) 

Zeff~  f  ~kRR^  Sin  ~  J  Jn(kR)  cos  «(<M0 
K  n=i  r  *  12  r„»  /_rtpM2 


t«tt>P-[*»M5> 


(d)  Computer  Program 


Using  Eqs.  (31)  and  (32),  Wu  and  Rosenbaum  have  solved  for  and  Ez  to  obtain  an 
input  wave  impedance  £2///<j>,  and  the  full  scattering  matrix  for  the  3-port  circulator: 


"lnwave 


=  -zd- 


±eff\j  D 

n  'lc?-c2c3 


Sn  =  522  =  S33  =  1  +  --^C)  ^C3  H 

jZeff{D )  j 

S21  =  S13  =  532  =  — •  (39) 

y  Zeff(  D ) 

S31  =S23=5,2= 

JZeff{D) 

Here  is  a  normalizing  impedance  (of  a  plane  wave  in  the  medium  just  outside  the  ferrite 
disk), 

zd=  2,11  Q  l(zd)m,  (40) 

where  Cj  is  the  relative  dielectric  constant  of  the  region  just  outside  the  disk,  and 

D=C}+  C\  +  C33  -  3  Ci  C2  C3  (41) 

The  values  of  Cl,  C2,  and  C3  are 


C1=3L^+  i  [sin VF j AgAs. -  JEZuL. 
2A0  n=l‘  n2xF  I  Dn  J2Zeff 


C2  =  —  B0  +  y  [sin^F 

2  2Ao  n 24/  , 

=  V  fanfoE 

3  2A0 


&n  COS 

Ft*) 

fi  flK 

tV*J 

sin 

ft) 

1  Dn 

^  cos 

[2  n  ffj 

♦ 

y_aK_ 

fl  2  sin 

1  Dn 
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where 


°n  =  A}  -  ,  with  An  =  Jn(kR)  and  Bn  =  Jn(kR) . 

This  infinite  series  solution  is  not  actually  correct  since  it  assumes  an  abrupt  step  of  the 
RF  fields  at  the  edges  of  the  three  input/output  transmission  lines.  The  abrupt  field  change 
cannot  actually  happen;  therefore,  some  number  of  terms  less  than  infinity  is  more  nearly 
accurate.  Bosma  [1]  speaks  of  a  simulation  using  18  terms,  Wu  and  Rosenbaum  [3]  use 
three  terms,  Schloemann  and  Blight  [4]  use  nine  terms.  In  many  test  calculations  using  the 
program  developed  here  it  has  been  observed  that  less  than  5%  difference  exists  in  the 
performance  parameters  between  using  three  terms  and  using  nine  terms.  The  Wu  and 
Rosenbaum  [3]  equations  can  be  programmed  almost  directly.  The  input  impedance, 
however,  should  be  denormalized  from  and  instead  normalized  to  the  characteristic 
impedance  of  the  input/output  stripline  or  microstrip  lines.  Also,  the  sign  of  the  input 
reactance  must  be  reversed  as  indicated  by  SaJay  and  Peppiatt  [6]. 

(e)  Computer  Program  Allowing  Negative  Effective  Permeability 

Most  circulators  above  about  2  GHz  are  biased  at  an  external  magnetic  field  in  the 
vicinity  of  4 jcMs,  which  is  a  field  level  below  the  one  required  to  produce  the  ferromagnetic 
loss  peak  at  or  near  the  circulator  center  frequency;  this  is  referred  to  as  below  resonance 
operation.  The  larger  field  necessary  to  operate  above  resonance  is  generally  not  physically 
practical  above  2  GHz.  The  internal  field  values  are  small  in  a  ferrite  circulator  disk  that  is 
biased  below  resonance  because  of  the  demagnetizing  effect  of  the  geometrical  shape  of  the 
disk.  The  demagnetizing  factor  with  the  bias  magnetic  field  perpendicular  to  a  disk  is 
AkMs  ,  so  that 


tfint  =  #ext  ~  AkMs  .  (43) 

With  relatively  small  internal  bias  field,  and  with  operating  frequency  far  above  the 
ferromagnetic  resonance  frequency,  the  imaginary  parts  of  \i  and  k  approach  zero,  the  real 

part  of  fi  approaches  unity,  and  the  real  part  of  Jcis  approximately  [-  (y/2  n)  •  4  n  Ms/f], 
as  shown  in  Fig.  5  of  Ref.  [2]. 

Schloemann  and  Blight  [4]  note  that  when  the  magnitude  of  v  becomes  larger  than 

unity,  which  happens  at  about  fm  =  (7/2tt  )  •  4 nMs,  then  fieff  =  (fi2  -  k2)/  fi  becomes 
negative.  To  make  calculations  possible  below  fm  they  rederived  the  Green’s  function 
terms  for  faff<  0.  The  result  is  that  V0,  Vn,  and  Un  of  the  new  Green’s  function,  written 
as 


G  (  0  ;  <j>')  =  V0  +  I  V„  cos  n(  <J>  -  <J>’)  +  Un  sin  n{  <j>  -  <|>’),  (44) 

n=l 

change  such  that  each  sign  is  reversed,  each  Bessel  function  Jn(kR)  becomes  a  modified 
Bessel  function  In(kR )  of  the  same  kind  and  order,  each  Bessel  function  derivative 
becomes  a  modified  Bessel  function  derivative  of  the  same  kind  and  order,  and  the  positive 
of  faff  is  used. 

To  keep  the  straightforward  computer  program  approach  used  by  Wu  and  Rosenbaum 
to  solve  the  circulator  problem,  it  is  necessary  to  find  relations  between  the  Green’s 


9 


R.  E.  NEIDERT 


function  terms  given  in  Refs.  [1]  and  [4]  and  the  Cl,  C 2,  and  C 3  terms  of  Ref.  [3]. 

Writing  Eq.  (37)  with  <J>  =  -  */3  and  <J>'  =  -  fl/3  corresponding  to  Port  1  of  the  circulator, 
with  r  =  R, 


But 


Therefore 


jjcZd 

2  Zeff' 


(45) 


(46) 


(47) 


This  expression  shows  how  Cl  is  modified  when  the  terms  of  G  (~nj 3  ;  -it] 3)  are  modified 
for  f<fm •  Similar  correlations  exist  between  C2  and  G  {-nj 3  ;  +7T/3)  ,  and  between  C3 
and  G  (-n/3  ;  it),  as  follows: 


and 

C3=/^[c(-f;K)].  (49) 

CIRCREN  uses  these  modified  expressions  in  its  calculations. 

There  is  a  small  difference  between  Refs.  [3]  and  [4]  in  which  Ref.  [3]  uses 

(sin2nTVn2  T^in  its  summation  of  terms  where  Ref.  [4]  uses  *F.  In  the  limit  the 
summation  term  approaches  'P  so  that  the  two  are  interchangeable  if  n  is  reasonably  large. 
Differences  between  the  two  were  very  minor  in  calculations  that  use  n  from  3  to  18. 

CIRCREN  is  set  up  to  use  the  (sin2  nx¥/n2'¥)  summation. 

3.  COMPUTER  PROGRAM  OPERATION 

The  theory  described  has  been  programmed  in  FORTRAN.  The  program  is  named 
CIRCREN.  Appendix  A  lists  the  February  1992  version  of  the  program  and  a  sample  run. 
It  is  currently  set  up  for  below  resonance  calculations,  for  applications  above  2  GHz 
assuming  that  the  bias  field  is  approximately  4 nMs.  It  calculates  the  characteristic 
impedance  of  stripline  or  microstrip  transmission  lines  in  and  out.  The  number  of  Green's 
function  terms  is  left  to  the  user  to  supply,  but  3  to  18  is  a  reasonable  range.  The  user 
supplies  the  4 nMs,  relative  dielectric  constant,  and  radius  of  the  circulator  puck,  and  the 
dielectric  constant  of  the  material  outside  the  ferrite  disk,  along  with  the  overall  height  and 
the  width  of  the  three  input/output  transmission  lines.  The  subfunction  BF  calculates  the 
Bessel  functions  directly  from  a  series,  switching  to  modified  Bessel  function  calculation 


10 


NRL  REPORT  9381 


when  necessary  by  means  of  the  KEY  value.  Recursion  relations  are  used  to  calculate  the 
Bessel  function  derivatives. 

Virtually  exact  agreement  (with  the  sign  of  the  input  reactance  reversed)  has  been 
acv  ‘eved  with  Wu  and  Rosenbaum's  calculations  by  using  n  =  3,  and  with  Schloemann  and 
Blight's  calculations  with  negative  jj.eff  and  with  n  =  9,  as  shown  in  Figs.  2  and  3, 
respectively.  Good  agreement  with  Salay  and  Peppiatt's  calculations,  compared  with  their 
measured  results,  is  shown  in  Fig.  4;  they  conclude  that  an  effective  radius  somewhat  less 
than  the  actual  radius  would  produce  better  calculated  agreement  with  the  measured  results. 
These  measured  data,  which  agree  with  their  own  independent  derivation,  are  also  their 
corroboration  that  the  sign  of  Bosnia's  reactive  term  is  reversed  from  the  correct  sign.  The 
calculations  of  CIRCREN  use  the  sign  of  the  reactance  that  agrees  with  Salay  and 
Peppiatt's  calculations  and  measurements.  The  low-frequency  extension  is  believed  to  be 
valid  only  down  to  about  f,J 2  as  discussed  in  Ref.  [4]. 


— O—  REACT  (WU  &  RO)  — O—  RES  (WU  &  RO) 

--O  REACT  (NEIDERT)  RES  (NEIDERT) 


Fig.  2  -  Comparison  with  Wu  and  Rosenbaum 


— O—  INSERTION  LOSS  (S  &  B)  — o-  ISOLATION  (S  &  B) 

--O-  INSERTION  LOSS  (NEIDERT)  *  -»  ISOLATION  (NEIDERT) 


R.  E.  NEIDERT 


Fig.  4  -  Comparison  with  Salay  and  Peppiatt 


4.  SUMMARY 

The  main  objective  here  has  been  to  produce  a  computer  program  to  calculate  the 
intrinsic  performance  of  a  Y-junction  circulator  from  physical  dimensions  and  magnetic 
properties  of  the  ferrite  material  and  embedding  structure.  This  has  been  done,  with  the 
program  set  up  for  the  most  common  case  of  below  resonance  operation.  The  mathematical 
manipulation  described  in  Section  2(e)  permitted  the  programming  of  the  convenient 
equation  forms  described  in  Section  2(d). 

Excellent  agreement  with  other  authors'  calculations  and  with  measured  data  in  a 
carefully  controlled  experiment  [6j  has  been  shown.  As  a  result,  we  can  confidently  say 
that  this  program  can  be  used  for  further  work  on  optimization  of  circulator  bandwidth, 
size,  and  monolithic  circuit  compatibility. 
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CIRCREN  CALCULATES  THE  PERFORMANCE  OF  A  FERRITE  CIRCULATOR 
INCLUDING  DIELECTRIC  LOSS, 

WITH  MU-EFFECTIVE  ALLOWED  TO  BE  NEGATIVE. 

COMPLEX  Cll, Cl , C2, C3,C4,C5,C6,C7 
COMPLEX  SM, AL1, AL2, AL, BE, GA, ZIN, ZINWU, YIN 
COMPLEX  EFC, ZEFF, SIC, ZIBCWU, ZIBC, SR, BO, AO 
COMPLEX  Z,  X9, Y9, RJ9, T9, S9, BF,  CIA,  C1C 
COMPLEX  BN, AN, CK2, CK22, C1BD, C1B,  C2BR,  C2BI ,  AL1 1 

C=2. 9979E10 
GAM=2 . 8E6 
PI=3. 141593 
U0=1 . 2566E-8 
E0=8 . 8543E-14 

50  WRITE (*,*) 'FOR  STRIPLINE  ENTER  1  ,  FOR  MICROSTRIP  ENTER  2 
READ ( * , * ) kstmi 

IF (KSTMI .EQ. 1)  GO  TO  26 
IF (KSTMI .EQ. 2)  GO  TO  26 
WRITE (*,*) 'SL/MS  CHOICE  NOT  1  OR  2' 

GO  TO  50 

26  WRITE (*,*)' ENTER  NUMBER  OF  GREENS  FUNCTION  TERMS  TO  USE 
READ (*, * ) NGFT 

WRITE (*,*) 'ENTER  4PIMS (G) , DIEL  CON  FERR, LOSS  TAN  FERR  :' 

READ (*,  *) FPMS,  EF, TAND 

WRITE (*,*) 'ENTER  DIEL  CON  OUTSIDE  FERR  :' 

READ (* , * ) ED 

FM=GAM*FPMS 
ZD=120 . *PI / SQRT (ED ) 

WRITE (*,*) 'ENTER  DIMS  RADIUS, LINE  W,DIEL  TH, LINE  TH  (ALL  CM) 
READ(*,*)R,W,B,T 

IF (KSTMI . EQ . 1 ) THEN 

C  CALCULATE  STRIPLINE  Z0 

IF (T . EQ . 0 . ) T=B/ 1000000 . 

X=T/B 

AM=2./(1.+2.*X/(3.*(1.-X))) 

DWBTA=X/PI/ (l.-X) 

DWBTB= (X/ (2. -X) ) **2+ ( ,0796*X/ ( (W/B) +1. 1*X) )  **AM 
DWBT=DWBTA* (1 . 5*LOG (DWBTB) ) 

WPBT= (W/ (B-T) ) +DWBT 

FSLA=  (8  .  /PI/WPBT)  +SQRT(  (8  .  /PI/WPBT)  **2  +  6.27) 

C  WU'S  STRIPLINE  GEOMETRIC  FACTOR, F, FOLLOWS 

FSL= ( . 5/PI) *LOG  fl .  +  (4 . /PI/WPBT) *FSLA) 

Z0SL= (ZD/2 . ) *FSL 
Z0=Z0SL 

ELSE 

C  CALCULATE  MICROSTRIP  Z0 

ETA0=120 . *PI 
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IF (T . EQ . 0 . )T=B/1.E6 

IF ( (W/B) .LE. (1./2./PI) ) THEN 

WPB= (W/B) + (1 . 25/PI) * (T/B) * (1 . +LOG ( 4 . *PI*W/T) ) 

ELSE 

WPB= (W/B) + (1.25/PI)* (T/B) * (1 . +LOG (2 . *B/T) ) 

ENDIF 

IF ( (W/B) .LE. 1 . ) THEN 

FWHL1= (l.+12.*B/W) ** (-. 5) +. 04* (1 .-W/B) **2. 

EE= ( (ED+1 . ) /2 . ) + ( (ED-1 . ) /2 . ) *FWHL1- (ED-1 . ) *T/B/4 . 6/SQRT (W/B) 
Z0MS=  (ETAO/  (2.  *PI*SQRT  (EE)  )  )  *LOG  (  (8  .  /WPB)  -i 
+  .-25*WPB) 

ELSE 

FWHG1= (1 . +12 . *B/W) ** (-. 5) 

EE=( (ED+1 . ) /2. ) + ( (ED-1 . ) /2. ) *FWHG1- (ED-1 . ) *T/B/4 . 6/SQRT (W/B) 
Z0MS= (ETAO/SQRT (EE) ) / (WPB+1 . 393+ . 667*LOG (WPB+ 

+  1.444)) 

ENDIF 

Z0=Z0MS 

ENDIF 

WRITE (*, *) ' Z0=  ' , ZO 
SI=ASIN(W/2./R) 

WRITE (*,*) 'ENTER  (IN  GHZ):  FLOW, FHIGH, FSTEP ' 

READ (* , * ) FREQL, FREQH,  FREQS 

WRITE (* , * ) '  FREQ  RIN  XIN  RTN  LS  ISOL  INS 

+LS  GIN  BIN' 


F=FREQL*1.E9 


200  KEY=2 

EFC=CMPLX (EF, -EF*TAND ) 

IF (ABS (F-FM)  .LT.  .001E9)  F=1.001*FM 

AK=-FM/F 

U=1 . 

UEFF= (U*U-AK*AK) /U 

IF (F  .LT.  . 999*FM) THEN 
UEFF=-UEFF 
KEY  =  1 

EFC=CONJG (EFC) 

ENDIF 

ZEFF=SQRT (U0*UEFF/E0/EFC) 

SIC1=PI/SQRT (3 . ) /I . 84 
SIC=SIC1*ZD*ABS (AK/U) /ZEFF 

ZIBCWU=2 . *ZD/ (1.+ (SIC/SI) * (SIC/SI) ) 
ZIBC=ZIBCWU*ZO/ZD 

SR=2 . *PI*F*SQRT (UEFF*EFC) *R/C 

B0=BF (0, SR, KEY) 

A0=-BF (1, SR, KEY) 

IF (KEY . EQ. 1 ) A0=-A0 


17 


R.E.NEIDERT 


ClA=SI*B0/2. /AO 
ClC=PI*ZD/2. /ZEFF 

AL11=CMPLX (0  .  ,  1.  ) 

Cll =C1A+AL11*C1C 
IF(F.LT.FM)  C11=-C1A+AL1 1 *C1C 


C1B=0 . 

C2BR=0 . 

C2BI=0. 

DO  100  N=l, NGFT 
RN=N 

BN=BF (N, SR, KEY) 

AN= (SR*BF (N-l, SR, KEY) -RN*BF <N, SR, KEY) ) /SR 

CK1= (SIN (RN*SI ) *SIN (RN*SI) ) / (RN*RN*SI) 

CK2=RN*AK/U/SR 

CK22=CK2*CK2 

C1BD=AN*AN-CK22*BN*BN 

C1B=C1B+CK1*AN*BN/C1BD 

CK3=COS (2. *RN*PI/3. ) 

CK4=SIN  <2 . *RN*PI / 3 . ) 

C2BR=C2BR+CK1 * (AN*BN*CK3) /C1BD 
C2BI=C2BI+CK1*CK2*BN*BN*CK4 /C1BD 
100  CONTINUE 

IF(F.LT.FM)  C1B=-C1B 
C1=C11+C1B 

C2=C1A+C2BR-AL1 1 *C2BI 
C3=C1A+C2BR+AL11*C2BI 
IF (F . LT . FM) THEN 
C2=-C2 
C3=-C3 
ENDIF 

C4=C1*C1-C2*C3 

C5=Cl*Cl*Cl+C2*C2*C2+C3*C3*C3-3 . *C1*C2*C3 

C6=C2*C2-C1*C3 

C7=C3*C3-C1*C2 

SM=PI*ZD/ZEFF/C5 

AL1=CMPLX (0 .  ,  -1 . ) 

AL2=SM*C4 
AL=1 . +AL1*AL2 
BE=AL1*SM*C6 
GA=AL1*SM*C7 

ZINWU=-ZD+AL1*2 . *ZEFF*CS/PI/C4 
ZIN=ZINWU*ZO/ZD 

RIN=REAL (ZIN) 

XIN=AIMAG ( ZIN) 

C  SALAY/PEPPIATT  REACTANCE  SIGN  CORRECTION 

XIN=-XIN 


DBRTNL=20 . *LUG10 (CABS (AL) ) 

DUISOL=20 . *LOG10 (CABS (BE) ) 

DBINSL=20 . *LOG10 (CABS (GA) ) 

YIN=1 . /ZIN 
GIN=REAL (YIN) 

BIN=AIMAG (YIN) 

C  SALAY/PEPPIATT  ADMITTANCE  SIGN  CORRECTION 

BIN=-BIN 
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WRITE (*,  10 ) F/l . E9, RIN, XIN, DBRTNL, DBISOL, DBINSL, GIN, BIN 
10  FORMAT (IX, 6 (IX, F8. 2) , 2(1X,F8.4) ) 

F=F+FREQS*1 .E9 

IF (F . GT . FREQH* 1.00001E9) GO  TO  210 
GO  TO  200 
210  STOP 

END 

COMPLEX  FUNCTION  BF(IP,Z,KEY) 

IF (IP . GT. 20 ) GOTO  1450 
IF (Z . GT . 20 . ) GOTO  1450 
X9=Z/2 . 

Y9=-X9*X9 

IF (KEY . EQ . 1 )  Y9=X9*X9 

J9=l 

RJ9=J9 

IF9=50 

IF (IP . EQ . 0 ) GOTO  1340 
DO  1000  K9=l , IP 
RK9=K9 

RJ9=RJ9*X9/RK9 
1000  CONTINUE 

1340  T9=RJ9 

DO  2000  K9=l , IF9 
RK9=K9 
RIP=IP 

RJ9=RJ9*Y9/RK9/ (RIP+RK9 ) 

S9=T9+RJ9 

IF (S9 . EQ. T9) GOTO  1470 
T9=S9 

2000  CONTINUE 

WRITE (*,*} 'MAX  NO  BESSEL  CALC  ITERATIONS  EXCEEDED.' 

WRITE (*,*)' SELECT  NEW  F, EF, 4PIMS, OR  R' 

STOP 

1450  WRITE (*,*) 'THE  VALUE  OF  BESSEL  ORDER  OR  ARGUMENT  GT  20.' 
WRITE (*,*) 'SELECT  NEW  F, EF, 4PIMS, OR  R.’ 

STOP 

1470  BF=T9 
RETURN 
END 
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Output  from  CIRCREN 

FOR  STRIPLINE  ENTER  1  ,  FOR  MICROSTRIP  ENTER  2  : 

1 

ENTER  NUMBER  OF  GREENS  FUNCTION  TERMS  TO  USE  : 

9 

ENTER  4PIMS (G) , DIEL  CON  FERR,  LOSS  TAN  FERR  : 

1780, 14,0 

ENTER  DIEL  CON  OUTSIDE  FERR  : 

10 

ENTER  DIMS  RADIUS, LINE  W,DIEL  TH, LINE  TH  (ALL  CM): 
.3, .417,  .125,0 
Z0=  7.862183 

ENTER  (IN  GHZ):  FLOW, FHIGH, FSTEP 
2,20,  .5 


FREQ 

RIN 

XIN 

RTN  LS 

I  SOL 

INSLS 

GIN 

2.00 

6.57 

0.90 

-19.27 

-18.26 

-0.12 

0.1493 

2.50 

6.79 

1.09 

-19.68 

-18.72 

-0.11 

0.1436 

3.00 

7.06 

1.23 

-20.17 

-19.29 

-0.09 

0.1374 

3.50 

7.39 

1.33 

-20.73 

-19.96 

-0.08 

0.1311 

4.00 

7.76 

1.35 

-21.29 

-20.70 

-0.07 

0.1250 

4.50 

8.17 

1.27 

-21.77 

-21.42 

-0.06 

0.1195 

5.00 

8.59 

1.09 

-22.02 

-22.01 

-0.05 

0.1146 

5.50 

8.98 

0.77 

-21.89 

-22.26 

-0.05 

0.1105 

6.00 

9.30 

0.32 

-21.34 

-22.00 

-0.06 

0.1074 

6.50 

9.49 

-0.24 

-20.44 

-21.23 

-0.07 

0.1053 

7.00 

9.52 

-0.87 

-19.34 

-20.14 

-0.09 

0.1041 

7.50 

9.37 

-1.51 

-18.17 

-18.90 

-0.12 

0.1041 

8.00 

9.02 

-2.10 

-17.01 

-17.65 

-0.16 

0.1051 

8.50 

8.52 

-2.58 

-15.88 

-16.47 

-0.22 

0.1075 

9.00 

7.91 

-2.93 

-14.78 

-15.37 

-0.28 

0.1112 

9.50 

7.22 

-3.12 

-13.68 

-14.33 

-0.36 

0.1166 

10.00 

6.48 

-3.24 

-12.41 

-13.31 

-0.48 

0.1235 

10.50 

5.41 

-3.64 

-9.92 

-11.97 

-0.79 

0.1272 

11.00 

4.93 

0.20 

-12.79 

-12.06 

-0.53 

0.2024 

11.50 

5.10 

-0.86 

-13.03 

-11.22 

-0.58 

0.1907 

0.2090 

12.00 

4.75 

-0.43 

-12.07 

-10.17 

-0.75 

12.50 

4.55 

0.36 

-11.42 

-9.01 

-0.96 

0.2184 

13.00 

4.58 

1.45 

-10.86 

-7.64 

-1.27 

0.1983 

13.50 

5.03 

2.94 

-10.21 

-6.04 

-1.83 

0.1483 

14.00 

6.44 

5.00 

-9.30 

-4 . 34 

-2.89 

0.0969 

14.50 

10.65 

7.29 

-8.13 

-2.85 

-4 . 86 

0.0640 

15.00 

19.91 

3.26 

-7.01 

-1.95 

-7.87 

0.0489 

15.50 

16.57 

-9.45 

-6.19 

-1.79 

-10.13 

0.0455 

16.00 

8.50 

-10.19 

-5.52 

-2.17 

-9.50 

0.0483 

16.50 

1.28 

-7.76 

-1.43 

-7.10 

-10.71 

0.0207 

17.00 

4.26 

-5.79 

-5.90 

-2.76 

-6.70 

0.0826 

17.50 

3.18 

-4.60 

-5.22 

-3.40 

-6.16 

0.1018 

18.00 

2.50 

-3.47 

-4.66 

-3.98 

-5.88 

0.1365 

18.50 

1.95 

-2.40 

-3.99 

-4.65 

-5.89 

0.2036 

19.00 

1.30 

-1.25 

-2.83 

-5.88 

-6.58 

0.4000 

19.50 

0.04 

1.30 

-0.09 

-19.96 

-20.19 

0.0238 

20.00 

4.47 

2.75 

-9.22 

-3.66 

-3  47 

0.1623 

STOP 


BIN 
0.0205 
0.0229 
0.0240 
0.0235 
0.0217 
0.0186 
0.0145 
0.0094 
0.0037 
0.0027 
0.0096 
0.0168 
0.0245 
0.0326 
0.0411 
0.0505 
0.0617 
0.0857 
0.0081 
0.0323 
0.0188 
0  0175 
0.0629 
0.0866 
0.0752 
0.0438 
0.0080 
0.0260 
0.0579 
0.1255 
0.1120 
0.1471 
0.1897 
0.2509 
0.3841 
0.7708 
0.1000 


20 


